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Abstract 

We give a variational proof of global stability for bistable travelling waves of 
scalar reaction-diffusion equations on the real line. In particular, we recover some 
of the classical results by P. Fife and J.B. McLeod (1977) without any use of the 
maximum principle. The method that is illustrated here in the simplest possi- 
ble setting has been successfully applied to more general parabolic or hyperbolic 
gradient-like systems. 



1 Introduction 

The purpose of this work is to revisit the stability theory for travelling waves of reaction- 
diffusion systems on the real line. We are mainly interested in global stability results 
which assert that, for a wide class of initial data with a specified behavior at infinity, the 
solutions approach for large times a travelling wave with nonzero velocity. In the case of 
scalar react ion- diffusion equations, such properties have been established by Kolmogorov, 
Petrovski & Piskunov [TT], by Kanel [9j [TO], and by Fife & McLeod [H [5] under various 
assumptions on the nonlinearity. The proofs of all these results use a priori estimates and 
comparison theorems based on the parabolic maximum principle. Therefore they cannot 
be extended to general reaction-diffusion systems nor to scalar equations of a different 
type, such as damped hyperbolic equations or higher-order parabolic equations, for which 
no maximum principle is available. However, these methods have been successfully applied 
to monotone reaction-diffusion systems [15], HE], as well as to scalar equations on infinite 
cylinders [TH [16] . 

Recently, a different approach to the global stability of bistable travelling waves has 
been developped by the second author [13J. The new method is of variational nature and 
is therefore restricted to systems which admit a gradient structure, but it does not make 
any use of the maximum principle and is therefore potentially applicable to a wide class 
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of problems. The goal of this paper is to explain how this method works in the simplest 
possible case, namely the scalar parabolic equation 



u t = u xx - F'(u) , (1) 

where u = u(x,t) G R, x G R, and t > 0. We shall thus recover the main result of 
Fife & McLeod pi] under slightly different assumptions on the nonlinearity F, with a 
completely different proof. The present article can also serve as an introduction to the 
more elaborate work [13] , where the method is developped in its full generality and applied 
to the important case of gradient reaction-diffusion systems of the form u t = u xx — W(u), 
with u G R n and V : R n —>■ R. A further application of our techniques is given in [7], 
where the global stability of travelling waves is established for the damped hyperbolic 
equation au tt +u t — u xx — F'(u), with a > 0. 

We thus consider the scalar parabolic equation (CQ), which models the propagation of 
fronts in chemical reactions [2], in combustion theory 0,(10], and in population dynamics 
[U E]. We suppose that the "potential" F : R — ► R is a smooth, coercive function with 
a unique global minimum and at least one additional local minimum. More precisely, we 
assume that F G C 2 (R) satisfies 

liminf uF'(u) > . (2) 

\u\— too 

In particular, F(u) — > +oo as |u| — > oo. We also assume that F reaches its global 
minimum at u — 1: 

F(l) = -A < , F'(l) = , F"(l) > , (3) 
and has in addition a local minimum at u = 0: 

F(0) = F'(0) = , F"(0) = p > . (4) 
Finally, we suppose that all the other critical values of F are positive, namely 

jw G R F'(u) = , F(u) < o} = {0 ; 1} . (5) 
A typical potential satisfying the above requirements is represented in Fig. 1. 

Under assumptions (jSJ)-©, it is well-known that Eq.(JT]) has a family of travelling waves 
of the form u(x,t) = h(x — c*t) connecting the stable equilibria u — 1 and u — 0. More 
precisely, there exists a unique speed C* > such that the boundary value problem 

h"(y) + cM(y) - F'(h(y)) = , y G R , 

/i(-oo) = 1 , /i(+oo) = , 1 > 

has a solution h : R — > (0, 1), in which case the profile h itself is unique up to a translation. 
Moreover h G C 3 (R), h'(y) < for all y G R, and converges exponentially to its limits 
as y — > ±oo. 

This family of travelling waves plays a major role in the dynamics of Eq.flT]), as is 
shown by the following global convergence result: 
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Fig. 1: The simplest example of a nonlinearity F satisfying assumptions ©-I©. 



Theorem 1.1 Let F G C 2 (M) satisfy assumptions (d))-([3p. Then there exist 5 > and 
v > such that, for all initial data Uq G C°(M) with 



limsup \uq(x) — 1| < 5 , limsup |u (x)| < 5 , 



(7) 



Eq. (QJ) has a unique global bounded solution satisfying u(x,0) = Uq(x) for all x G M. In 
addition, there exists x G M such that 



sup 



u(x, t) - h(x - c*t - x ) = 0(e~ ut ) 



as t 



-oo . 



Theorem 11.11 was first proved by Fife & McLeod [HE] under the additional assumption 
that < u (x) < 1 for all In that case u(x, t) G [0, 1] for all x G M and all t > 

by the maximum principle, so that the coercivity assumption ([2]) is not needed. As is 
mentioned in [3] , the results of [4] can be extended to arbitrary initial data satisfying ([7]) 
provided that uF'(u) > for all u ^ [0, 1], a condition that is more restrictive than ([2]) 
in the sense that F is not allowed to have critical points outside the interval [0, 1]. The 
simplest case considered in jl] is when F has exactly one critical point in the open interval 
(0,1), a situation in which condition ([5]) is clearly met. However, Fife & McLeod also 
study the case where F has three critical points in the open interval, including a local 
minimum at u = u* G (0, 1). In this situation there exists a travelling wave solution of 
([T]) with speed C\ > connecting u = 1 to u = u*, and also a travelling wave with speed 
c 2 G M connecting u = to u = 0. If C\ > c 2 , which is always the case if §5§ holds, 
there exists c* G (0, 1) such that §6§ has a solution h : M — > (0, 1), and the conclusion of 
Theorem 11.11 is still valid. If c\ < c 2 , there exists no travelling wave connecting u = 1 to 
u = 0, and the solution of (CQ) with initial data satisfying ([7|) converges as t — > oo to a 
superposition of two travelling waves [1]. 

Theorem 11.11 is a particular case of the general results obtained in [13], see Theorem 4 
in Section 9.6 of that reference. Therefore, there is no need to give here a complete proof. 
Instead we shall prove the convergence result ([HD under the additional assumption that 
the initial data u (x) decay rapidly to zero as x — » +oo. It is intuitively clear that the 
precise behavior of uq(x) near x = +oo should not play an important role, because the 
equilibrium u = ahead of the front is stable (this is in sharp contrast with the case of a 



3 



monostable front invading an unstable equilibrium, where the behavior ahead of the front 
is of crucial importance). However, this restriction allows to shortcut many technicalities 
and to give a much simpler proof in which the essence of the argument can be easily 
understood. 

Our approach is based on the fact that Eq.flJJ) possesses (at least formally) a gradient 
structure, not only in the laboratory frame but also in any frame moving to the right with 
a positive velocity. To see this, we introduce the following notation. If u(x, t) is a solution 
of flTJ, we define for any c > 

v(y,t) = u(y + ct,t) , or equivalent ly u(x,t) = v(x — ct,t) . (9) 

Setting y = x — ct we see that the new function v (y, t) satisfies 

Vf = Vyy + CVy ~ F'(v) . (lO) 

We now introduce the energy functional 

m = Je»(±v}+F(vj)dy, (11) 
and the corresponding energy dissipation functional 

V c [v}= ( e cy (v yy + cv y -F'(v)) 2 &y . (12) 

We also denote by H^{R) the Banach space 

Hi (R) = jt> G L°°(R) e cy/2 v G tf^R)} , (13) 

equipped with the norm \\v\\h^ = \\v\\l°° + \e cy l 2 v ||#i. Note that any v G Hl(M) decays 
to zero faster than e~ cy ^ 2 as y — > +oo. Since F(v) ~ f3v 2 /2 as v — >• by (BJ, it follows 
that £ c [v] < oo for all v G if^(R). Conversely, any v G L°°(R) such that v(y) — > as 
y — > +oo belongs to H^iM.) as soon as £^ c [f] < oo. 

If f (?/,£) is a solution of (JTUJ) with initial data v G ^(R), then v(-,t) G i^QR) for all 
t > and a direct calculation shows that 

^£M;t)] = -T>M;t)] < , t > . (14) 

In other words, the energy £ c is a Lyapunov function of system (1101) in if*(R). This 
observation is of course not new: in their original proof, Fife & McLeod [4] already used 
a suitable truncation of the functional £ c for the particular value c = c* to show that the 
solution v(y,t) of ffTUj) approaches a travelling wave for a sequence of times. However, 
the fact that Eq.(JTJ) has a whole family of (nonequivalent) Lyapunov functions has not 
been fully exploited until recently. The only reference we know where the implications 
of this rich Lyapunov structure are really discussed is a recent paper by Muratov [T2] . 
which contains a lot of interesting observations and a few general results concerning a 
wider class of systems than Eq.flT]), but fails to prove the convergence to travelling waves. 
The goal of the present article is to show that, in the simple case of Eq.flT]), the gradient 



4 



structure alone is sufficient to establish convergence, at least if we restrict ourselves to 
solutions which decay to zero rapidly enough as x — > +00 so that the energy functionals 
are properly defined. 

The main difficulty of this purely variational approach is that we do not have good a 
priori estimates on the solution v(y, t) = u(y + ct, t) in a moving frame with speed c > 0. 
First of all, it is not clear a priori that the energy S c [v(-,t)] is bounded from below (this 
will not be the case typically if c is too small) , and without this information it is difficult 
to really exploit the dissipation relation (fl4|) . Next, if we have a lower bound on S c [v(-, £)], 
we can deduce from ( fl4l) that the solution v(y,t) converges uniformly on compact sets, 
at least for a sequence of times, towards a stationary solution of ( TTUl) . but we cannot 
exclude a priori that this limit is just the trivial equilibrium v = (this will be the case 
typically if c is too large). To overcome these difficulties, the main idea is to track the 
position of the front interface in the following way. We fix positive constants /3 2 such 
that Pi < F"(0) < /3 2 , and we choose e > small enough so that 

A < F'\u) < p 2 , for all u G [-2e, 2e) . (15) 

Given a continuous solution of ([1]) satisfying the boundary conditions 

lim u(x,t) = 1 , lim u(x,t) = , t > , (16) 

X — > — 00 x— > + oo 

we define the invasion point x(t) as the first point starting from the right where the 
solution u(x,t) leaves an e-neighborhood of the equilibrium u = 0: 

x(t) = max jx G R \u(x, t)\ > ej . (17) 

In view of (jTBT) . it is clear that —00 < x{t) < 00 for all t > 0, and that \u(x(t),t)\ = e. A 
quantity similar to x(t) was also introduced in [12], where it is called the "leading edge". 

The strategy of the proof is to show that the solution u(x, t) converges uniformly on 
compact sets around the invasion point x(t) towards a suitable translate of the travelling 
wave ([6]). Using only the gradient structure, we can prove the following result: 



Proposition 1.2 Let F e C 2 (1R) satisfy assumptions (d))-([3]]. If uq G -f^c(IR) for some 
sufficiently large c > and m — 1 G i/ 1 (M_) ; then the solution u(x,t) of Eq.ljj\) with 
initial data Uq satisfies, for all L > 0, 

sup \u(x(t) + z,t) - h £ (z)\ >0, (18) 

ze[-L,+oo) 

where x(t) is the invasion point [F?\ ) and h £ is the travelling wave (0|) normalized so that 

h £ (0) = e. Moreover the map t 1— > x(t) is C l for t sufficiently large and x'(t) — > c* as 
t — > 00. 



As is explained above, the assumption u G H^W) is needed in order to use the energy 
functional S c without truncating the unbounded exponential factor e cy . The proof will 
show that is it sufficient to take here c > y/2A/e, where A is defined in (|3]) and e in (JTB1) . 
On the other hand, the assumption uq — 1 G if^R-) is just a convenient way to guarantee 
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that the first condition in (|T6|) is satisfied, but with minor modifications we can treat the 
more general case where \uq(x) — 1| is assumed to be small for large x < 0, as in ([7j). 

The local convergence established in Proposition 11.21 is the key step in proof of The- 
orem [LTJ Once (fT8j) is known, it remains to show that the solution u(x,t) converges 
uniformly to 1 in the region far behind the invasion point x(t). Such a "repair" is cer- 
tainly expected because u — 1 is the point where the potential F reaches its global 
minimum. A convenient way to prove this is to use a truncated version of the functional 



where F(u) = F(u) — F(l) > 0. In this way, we can show that the solution u(x,t) 
approaches uniformly onKa travelling wave (at least for a sequence of times), and using 
in addition the local stability results established in [T7] we obtain (JHD- We thus have: 

Corollary 1.3 Under the assumptions of Proposition [TTR there exist xq G R and v > 
such that (0|) holds. 

We conclude this introduction with a few comments on the scope of our method. 
First, it is clear that the assumptions (EJ)— (ISJ) are not the weakest ones under which 
Proposition ll.2l holds. A careful examination of the proof reveals that the only hypotheses 
that we really use are: 

HI: For all bounded initial data uo, Eq.^ has a (unique) global bounded solution. This 
is certainly true if ([2]) holds, but it is sufficient to assume, for instance, that F(u) — > 
+oo as \u\ — > oo, or that uF'(u) > whenever \u\ is sufficiently large. 

H2: F(0) = F'(0) = ; and there exists e > such that F"{u) > for all u G [-e,e\. 
This is automatically true if (jl]) holds, but u = need not be a strict local minimum 
of F. In particular Proposition 11.21 holds for the nonlinearities of combustion type 
considered in [9j [10] . 

H3: There exists a unique c > such that the differential equation v yy + cv y — F'{y) = 
has a bounded solution satisfying \v(0)\ = e, \v(y)\ < e for all y > 0, and v(y) — > 
as y — > +oo; furthermore, this solution is unique. Under assumptions (Jl])-©, we 
have c = c* and v = h e . In general, we can assume without loss of generality that 
that v is positive and converges to 1 as y — > — oo, so that F(l) < and F'(l) = 0. 
It also follows that F(u) > for all u < and that F has no critical point < 1 
with < 0. 

On the other hand, to prove that the solution of (CQ) given by Proposition 11.21 converges 
uniformly on K to a travelling wave we need the additional assumption: 

H4: There exists e' > such that the only bounded solution of the differential equation 
u xx — F'(u) = with \u(0) — 1| < e' is u = 1. This requires that F attains its global 
minimum at u — 1, and nowhere else. 





(19) 
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Finally, if we want the convergence to be exponential in time as in (151) , we need to assume 
that F"(l) > 0. 

Another comment concerns the variational structure of Eq.([T)). Due to the exponential 
weight e cy , it is clear that the energy functional £ c is not translation invariant. In fact, for 
any v G H^(M.) and any I G 1R, we have the relation £ c [v (■—£)] — e ce £ c [v}. This implies 
that the infimum of S c [v] is either or — oo. Under our assumptions on F, the transition 
between both regimes occurs precisely at the critical speed c* for which travelling waves 
exist: 

if c > c* , 



inf £ c \v] = i ... 

vdH\ (R) [ -OO it C < C 

Indeed, as was observed by Muratov [12] , for any c < c* + a/ c * + 4F"(0) we have the 
identity 

cS c [h] = (c-c*) / e^h'iyfdy, 

where h is the solution of ([6]). This shows in particular that S c [h] < when c < c*, hence 
inf £ c = — oo in that case. The fact that £ c > when c > c* is not obvious a priori, and 
will be established in the course of the proof of Proposition 11.21 see Corollary 14.31 Note 
also that £ c *[^] = 0, so that inf £ Ct = mm£ Ct = 0. 

The rest of the paper is organized as follows. In Section [21 we establish the basic 
inequalities relating the energy £ c , the dissipation V c , and the invasion point. Using these 
relations, we prove in Section [3] that the average speed of the invasion point x(t) has a 
limit Coo > as t — > oo. The core of the paper is Section @J where we show that Coo = c* 
and prove Proposition 11.21 The proof of Corollary 11.31 is then performed in the final 
Section [51 

Acknowledgements. The authors are indebted to S. Heinze, R. Joly, and C.B. Muratov 
for fruitful discussions. 



2 Preliminary estimates 

As the potential F is smooth and coercive, it is well-known that the Cauchy problem for 
the semilinear equation ^ is globally well-posed in the space of bounded functions, see 
e.g. [8J. Due to parabolic regularization, the solutions are smooth for t > and satisfy 
([I]) in the classical sense. Under assumption ([2]), one can also show that our system has 
a bounded absorbing set in the following sense: 

Lemma 2.1 There exists a constant B > depending only on F such that, for all initial 
data Uo G L°°(R) ; the (unique) solution u(x,t) of (TJP satisfies, for all sufficiently large 
t > 0, 

sup(\u(x,t)\ + \u x (x,t)\ + \u xx (x,t)\) <B. (20) 
Moreover, u(-,t) is bounded in Hf oc (M) for some s > 5/2 and all t > 1. 
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The uniform bound on \u(x, t) \ follows easily from the maximum principle, but it can 
also be established using localized energy estimates, see [T3J Section 9.1]. The bounds on 
the derivatives are then obtained in a standard way using parabolic regularization. 

From now on, we suppose that u$ G H^ o (R) for some Cq > (which will be specified 
later) and that u — 1 G Then the solution of with initial data Uq satisfies 

u(-,t) G H^(M) and u(-,t) - 1 G H l (R_) for all t > 0, because u = and u = 1 are 
(stable) equilibria of ([I]). In particular, the boundary conditions (jT6l) hold for all times, 
so that one can define the invasion point x(t) by (|T7|) . Also, since we are interested in the 
long-time behavior of u(x, t), we can assume without loss of generality that estimate ( |20l) 
is valid for all t > 0. 

As is explained in the introduction, we shall use the energy functionals £ c (for various 
values of c > 0) to prove that the solution u(x, t) converges to a travelling wave h locally 
around the invasion point x(t). A technical problem we shall encounter is that the invasion 
point, as defined in ( 1TT1) . need not be a continuous function of time and can therefore jump 
back and forth in an uncontrolled way. It is possible to avoid this difficulty using a more 
clever definition than (|T7|) . see [T5] , but we follow here another approach and just introduce 
a second invasion point defined by 



X(t) = maxjx G K \u(x,t)\ > 2e| . (21) 



Clearly, — oo < X(t) < x(t) < +oo for all t > 0. The important point is that an 
information on x at a given time provides an upper bound on X at later times: 

Lemma 2.2 There exists To > and Co > such that, for all to > 0, one has 

X{t) < x{t ) + C for all t G [t , t + T ] . (22) 

Proof: Fix to > 0. The solution of ([T|) satisfies 

u{t) = S(t - t )u(t ) - f S(t- s)F'(u(s)) ds = u x {t) + u 2 {t) , t > t , 

where S(t) = e td * is the heat semigroup. Take K > such that |F'(m)| < K whenever 
\u\ < B, where B is as in (120]) . Then ||u 2 (t)||L°° < K{t — to). On the other hand, by 
definition of x, we have \u(x, t )\ < e if x > x(t ) and \u(x, t )\ < B if x < x(t ). Using 
the explicit form of the heat kernel, we deduce that 

/ m 1 f (x ~ y)2 1 / \i , B / x — x(t ) \ 

ut(x,t) < / e k^ u(y,t ) dy < e + - erfc lM =) , 

J4n It - to) Jr 2 V /4 (t-to V 



'\<K{t - to) Jr " ' 2 \J4(t-t y 

where erfc(x) = (2/y/n) e~ z2 dz. We first choose T > such that KT < e/2, and 
then C > such that B erfc(C /\/4To) < e. Then, for all t G [to^o + ^o] an d all 
x > x(t ) + Co we have |u(a;,t)| < 2e, which implies that X(t) < x(t ) + C . □ 



We now derive the basic estimates on the energy (1TT1) and the energy dissipation ( 1121) 
which will be used throughout the proof. Given c G (0, Co), we define v (y, t) = u(y + ct, t) 
as in IQ, and we set 

E c (t) = £ c [v(-,t)} , D c (t) = V c [v(;t)} , t>0. (23) 
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Of course, v(y, t) depends also on the speed parameter c, but to simplify the notations this 
dependence will not be indicated explicitly. We also denote by y c (t), Y c (t) the invasion 
points in the moving frame: 

y c (t) = x(t) - ct , Y c (t) = X(t) - ct . (24) 

By construction, \v(y c (t), t)\ = e, \v(y, t)\ < e for y > y c {t) and \v(y, t)\ <2e for y > Y c {t). 
Remark that, by ( fl5l) . the following inequalities hold whenever \v\ < 2e: 



^v 2 < F{v) <^v\ fav 2 <vF\v) < p 2 v\ ft < F"(v) < fa . (25) 

Lower bound on E c : Using (fTTl) . (1251) . and the fact that F(u) > —A for all u E R, we 
find 



E c (t) = J e c «(-v 2 y (y,t)+F(v(y,t)))dy + j ^ (-v 2 y (y,t) + F(v(y,t))) dy 

rVc(t) poo -I n 

> / e c y(-A)dy+ M v 2 (y, t ) + ^v 2 (y,t))dy. (26) 

J-oc Jy c (t) KZ I J 

To estimate the last integral in (J26l) we recall that v(y c (t),t) 2 = e 2 , so that 

e cUt) e 2 = - T d y (e c yv 2 { yi t))dy = - f° (2v(y, t)v y (y, t) + cv 2 (y, t)) dy . 



Given d > — c, we have \2vv y \ < (c+d)v 2 + (c+d)~ 1 v 2 , hence 

e*<V< e c y(^-v 2 y (y,t) + dv 2 (y,t))dy . (27) 
Jut) Vc + « y 

If we choose d such that d(c + d) = f3\ and insert the resulting inequality into fl2"B"j) . we 
obtain 

£c(*) > e c5cW f-- + ke 2 ) , where < « < -(c + ^c 2 + 4ft) . (28) 
V c / 4 

This estimate shows in particular that the energy E c (t) is bounded from below as long 
as the invasion point y c (f) is bounded from above. Moreover, the lower bound is close to 
zero if y c (t) is large and negative. 

Variation of D c : It follows from ([TDD and (PJ that £> c (t) = / R e cy vf{y, t) dy. Differen- 
tiating this relation with respect to t and integrating by parts, we find 

~ D' c (t) = [ e<*(v t v tt )(y, t) dy = [ e cy v t {v tyy + cv ty - F"(v)v t ) dy 
1 Jr Jr 

= - [ e c yv 2 y (y,t)dy - [ e c y F"(v(y,t))v 2 (y,t) dy . 



Take C\ > such that F"(u) > —C±/2 whenever \u\ < B, where B is as in ( 120|) . Then 
the above relation shows that 

D' c (t) < C x D c {t) , t>0. (29) 
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This differential inequality implies that, if D c G L 1 (IR + ), then D c (t) — > as t — > oo. Since 
E' c (t) = —D c {t) by (fT4"j) . this will be the case as soon as E c is bounded from below. 

Lower bound on D c : Using (fl~2l again and integrating by parts, we find 

D c (t) = [ e cy (v yy + cv v -F'(v)) 2 (y,t)dy 
Jr 

= [ (v 2 yy + 2F"{v)v 2 y + F'(v) 2 ) (y, t) dy . 
Jr v ' 

We split the integration domain into (— oo, Y c (t)) and (Y c (t), +00). Using (|20|) . (1251) and 
the lower bound F"(v) > —Ci/2, we obtain 

A:(f) > -— e cY ^+ I e^(vl y (y,t) + 2p 1 vl(y,t)+py(y,t))dy. 



c 



Yc(t) 



Observe that, for any y G R, 

/oo 2 /*oo /-co 2 /-co 

e c ^ 2 dy > - e cy v 2 dy , and / dy > - e cy v 2 y dy . (30) 

„o -Ato Jya Jyo 

Indeed the first inequality is just fl27|) with e — 0, d — — c/2, and y c (t) replaced by ?/o> 
and the second inequality is similar. Thus, for any d > we have 

D c {t) > -5^!e*(*) + jT e^{(2/5 1 + ^-rf)^+(/? 1 2 + ^)t; 2 }( 2/ ,t)dy. (31) 



In an analogous way we find 



= / e c ^(-^( ? /,t) + F(t;( 2/ ,t))Jd Z /+ / e^(-^(y,t) + F( U (y,0) 



dy 



< -e*W+ r K^^ + V^ /))<!/;. (32) 



where X = (B 2 /2) + IT and A"' = sup{F(w) | \u\ < B}. If we now combine (131]), (j32j| and 
choose the particular value d = (3 2 — (P2 ~ Pi) 2 + c 2 /4) > 0, we arrive at 

D c (t) > 7 £ c (t)-^e c ^ , t>0, (33) 
c 

where 

< 7 < ~ (c2 2 t 4 zf/ » and ^2 = ^ 2 + 7^. 
2 + 4^2 

Inequality (1331) means that, if the invasion point Y c (t) is large and negative, the energy dis- 
sipation D c = —E' c is essentially proportional to the energy itself. This gives a differential 
inequality for E c (t) which, in view of Lemma [2 .2\ can be integrated as follows: 

E c (t) < e-^~^E c (t ) + ^ e c ^)+c ) ? t G ^ k + Tq] (34) 

c 

Remark 2.3 The constants Cq,C\ and Tq introduced in this section depend only on the 
potential F. In particular, they are independent of the solution u(x,t) and of the speed 
parameter c. Similarly, if we choose k = \ff3~i/2 and 7 = 2(31/(32, then the constants k, 7 
and C 2 depend only on F . 
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3 Existence of the invasion speed 



As in the previous section, we suppose that u(x, t) is a solution of (CQ) with initial data uo 
satisfying the assumptions of Proposition 11.21 We also assume that the bounds (1201) hold 
for all t > 0. If x(t) denotes the invasion point (ED, we define 

T • f 5 W I" ^ /or \ 

c_ = hmmt , c + = limsup . (ob) 

t^oo t ' t^oo t 

Our first result shows that the solution u(x,t) invades the stable equilibrium u = at a 
positive, but finite, speed. 

Proposition 3.1 One has c_ > and c + < oo. 

Proof: The proof relies on the lower bound ( 1281) . Assume that the initial data uq belong 
to Hl o (R) for some c > V2A/e, where A = -F(l) and e is as in flT5]). Using (TJED with 
c = c and k = c/2, we find that E c (t) > a e c5c(t) for some a > 0. Since £? c (t) < i? c (0) for 
alH > 0, it follows that y c (t) = x(t) — ct is bounded from above, hence 

Vc(t) 

c + — c + hm sup < c < oo . 

t^oo t 

On the other hand, since uq — 1 G and .F(l) = —A < 0, it is easy to verify 

that -Ec(O) = ^ c [mo] ~ —A/c as c — > 0. Thus if we take c > sufficiently small so that 
E c (0) < 0, it follows from ([28]) that > E c (0) > E c {t) > (-A/c) e c5c(i) for all t > 0. This 
implies that y c (t) = x(t) — ct is bounded from below, hence 

c_ = c + lim inf > c > . 

t— »oo t 

This concludes the proof. □ 
We next prove that the average invasion speed x(t)/t converges to a limit as t — > oo. 

Proposition 3.2 One /ias c_ = c+. 

Proof: We argue by contradiction. Assume that c_ < c+, and choose time sequences 
{t n } neN , {^}nGN such that t n -> oo, ^ -> oo and 

^(*n) g(*n) 

— ; — c - j > c + • 

t n^oo t n ra^oo 



Due to Lemma EH], upon extracting a subsequence we can assume that u(x(t n ) + z,t n ) 
converges in Hi oc (M) to some limit Woo(z). More precisely, for any L > 0, 

u(x(t n ) + z,t n ) > Woo(z) in H 2 ([-L,L}) , 

u t {x(t n ) + z,t n ) > WocW in L 2 ([—L,L]) , 



n— >oo 
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where w M G H^.(R) n L°°(R) and G satisfy Woe = w'^ - F'(woo). Moreover, 

by definition of the invasion point, we have |u>oo(0)| = e. 

Now, we fix any c G (c_, c+) and we observe that the invasion point y c (t) = x(t) — ct 
satisfies y c {t' n ) — * —oo and y c {t n ) — * +00 as n — > 00. Using first the lower bound fl28|) . we 
find 

# c (0 > — e c ^*«) ► , (36) 

C n— >oo 

hence (since E c is non-increasing) -E c (t) > for all t > 0. As ^(i) = —D c (t), we deduce 
that 

/>oo 

(37) 



/■oo 

/ D c (t) dt < E e (0) < 00 



and using in addition (I2"9"j) we conclude that -D c (t) — > as t — > 00. 
Next, we observe that, for all n G N, 

D c (t n ) = [ e cy v 2 (y,t n )dy = e c ^ [ e cz v 2 (y c (t n ) + z, t n ) dz 
Jr Jr 

= e c ^») / e cz (u t + cu x ) 2 (x(t n ) + z, t n ) dz . (3* 



Since D c (t n ) — > as n — > 00, the last integral in ( |38|) converges to zero as n — > 00, hence 
the limits u>oo,^oc satisfy Woo + cw'^ = 0. Incidentally, this means that w'^ + cw'^ — 
F'(woo) = 0, i.e. Woo is a travelling wave solution of (0Q) with speed c. Now the crucial 
point is that c G (c_,c+) is arbitrary. Obviously, the relation Woo + cw'^ = can be 
satisfied for two different values of c only if w'^ = 0, i.e. if Woo is identically constant. 
But then we must have F'faao) = 0, which is impossible in view of (125|) since \woo\ = e. 
This contradicts the assumption c_ < c + and concludes the proof. □ 



Remark 3.3 Another way to obtain a contradiction in the proof of Proposition ^. £1 which 
works even if u = is not a strict local minimum of F (see hypothesis H2 in the intro- 
duction), is to observe that the limiting function Woo{z) converges to zero as z — > +00. 
Indeed, proceeding as in l[2b)) . ( fffgj) and using [3ty) we find for all n G N: 



E c (t n ) = e c ^ f e cz ^-u 2 x + F(u))(x(t n )+ Zl t n )dz 



> e cy c (t„) y°° e cz (I u 2 + /| u 2 j ( - (tn) + ^ tn) dz--^ e c ^ ] (39) 

* JO c 

where k' = min{l, (c 2 +4/5i)(c 2 +4)" 1 }. As E c (t n ) < E c (0) and y c {t n ) — > +00 as n — > 00 ; 
we /iawe fry Fatou's lemma: 

1 

e c >^) 2 + «a^) 2 )^ 

/■°° 

< liminf / e C2 (w 2 + u 2 )(x(t n ) + z, t n ) dz < — < 00 . 

n-*oo J CK 

Thus Woo £ H^{R), and in particular w^z) — > as z — > +00. 27ms zs clearly impossible 
w 'oo = an ^ koo(0)| = e. 
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4 Local convergence to a travelling wave 



This section is devoted to the proof of Proposition [L2] Using the same notations as in the 
previous sections, we first prove that the solution u(x, t) of (JTJ converges for a sequence 
of times towards a travelling wave, locally in space around the invasion point. On this 
occasion we identify the invasion speed given by Proposition 13 . 21 with the unique speed c* 
for which travelling waves exist. 

Proposition 4.1 Let Cqo = c_ = c+. There exists a sequence t n — > oo such that, for all 
L > 0, 

/ e c °° z (u t + CooU x ) 2 (x(t n ) + z, t n ) dz > . (40) 

J-L n_>0 ° 

Proof: Since the left-hand side of ( 1401) is a nondecreasing function of L, it is sufficient 
to prove that, for any L > 0, there exists a sequence t n — > oo such that ( 140]) holds. We 
argue by contradiction and assume that there exist L > and <5 > such that 

J e c °" z (u t + c 00 u x ) 2 (x(t)+z,t)dz > 5 , (41) 

for all sufficiently large t. In fact, upon changing the origin of time, we can assume that 
(jjrj holds for all t > 0. In analogy with (T2~4"|) . we denote ?/(£) = x(t) — c^t. Two situations 
may occur: 




y(t) 



Fig. 2: If there exists a sequence t n — > oo such that is bounded from below, a contradiction is 

obtained by considering the dissipation of the energy E c in a moving frame with speed c > (c close 
to Cqq). If y(tn) > 1 for all n £ N, the set S c consisting of all n such that y(t n ) > (c — Coo)t„ increases as 
c — > Coo, and card(S' c ) — > oo. 



Case 1: There exists a time sequence t n — > oo such that {?/(£«) } n eN is bounded from 
below. Without loss of generality we can assume that t n+ i > t n + 1 and 7/(t n ) > 1 for all 
nGN (the second condition is easily achieved by translating the origin). 

Let K > be such that E c (0) < K for all c G [cqo, c ], where c > is as in the proof 
of Proposition 13.11 Take c G (c^, c ) sufficiently close to so that 

e ( c - Coo )L < 2 ) and ( c _ Coo )2 5 2 / eCoo2 dz < - , (42) 



4 
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where B is as in ([2D]) . Let y c (t) = x(t) — ct = y{t) — (c— c oc )t. Since > 1 for all 

n 6 N, it is clear that the cardinality of the set 

S c = {n G N | y c (tn) > 0} = {n G N | y(t n ) > (c - Coo )U 

becomes arbitrarily large Fig. 2. On the other hand, y c {t n ) — ► —oo as 

n — ► oo and this implies (as in the proof of Proposition 13. 21) that E c (t) > for all t > 0. 
But for all n G 5 C , we have by (1381). (T4"Tj), (|42l) 

A=(*n) = e c ^ [ e cz (u t + cu x ) 2 (x(t n ) + z,t n )dz 

/L x 
e cz (u t + cu x ) 2 (x(t n ) + z, t n ) dz > - , 

because e cz > \e CooZ for \z\ < L and (u t + cu x ) 2 > \{u t + c^Ux) 2 — (c — Coo) 2 !! 2 . Moreover, 
it follows from §M) that D c (t) > Djt n ) e~ Cl for all t G [t n - l,t„], hence 

E c (t n -1) - E c (t n ) = [ D c (t)dt > 5 -e~ c ' , n G S c . 

Jt n -l O 

If we choose c close enough to so that card(S' c ) > 8Ke Cl /5, we obtain a contradiction 
with the fact that E c (t) is positive, nonincreasing, and E C (Q) < K. 




Fig. 3: If y(t) — ► — oo a contradiction is obtained by considering the dissipation of the energy E Cn in a 
moving frame with speed c„ < on the time interval [0,t n ], where y(t„) = (c n — Coo)i n . We choose 
T 1 and then n large enough so that (cqq — c n )T < 1. 



Case 2: — ► — oo as t —>■ oo. In this case, there exists a sequence t n — ► oo such that 

y(t n ) < inf y(s) + 1 , for all n G N , (43) 

0<s<t n 

see Fig. 3. Indeed the function fi(t) = mf{y(s) | < s < t} is nonincreasing and /i(t) — > 
— oo as £ — > oo. For each n G N, we choose t n G [0, n] such that |/(t n ) < M 71 ) + 1- Then 
/x(n) < £t(i n )> hence (T431) holds. 
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Given some (large) n G N, we take c n < such that (c n — Coo)t„ = ?/(£„), or 
equivalently y Cn (tn) = 0. Since c n — > as n — > oo, we can assume that c n > Cqo/2 and 

e c °° 2 dz < - . (44) 
l 4 

If t G [0,t„], we have by (|13|) 

f7c(*) = + ( c oo - C n )t > y(t n ) - 1 + (Coo - c n )t = (c - Coo)(t n - t) - 1 . 

Using (|4T1) . (jHJ) and proceeding as in the previous case, we obtain 

D Cn (t) = ^ cn{t) [ e c " z (u t + c n u x ) 2 (x(t)+z,t)dz 



r 

> p Coc((c n -Coo)(tn-t)-l) _ 

8 ' 

hence for all T <t n : 

rt„ r 

/ A*(*)df > Te c "« Cn - Coo > T - 1 >~ . (45) 

On the other hand, there exists K > such that E Cn (0) < K for all n, and since y Cn (t n ) = 
we know from (1251) that E Cn (t n ) > —A/c n . Thus 

Ajt) dt = £ Cn (0) - £ c „(t n ) <ir + -<K + _. (46) 

Cn Coo 

If we now choose T > large enough so that Te~ 2c °°5 > 8K + WAfc^,, and then n£N 
large enough so that t n > T and (c^ — c n )T < 1, we obtain the desired contradiction by 
comparing (|45p and (I46p . □ 



Corollary 4.2 One /ias c_ = c + = c* ; and i/iere exists a sequence t n — > oo snc/i t/iat, /or 
a// L > 0, 

SUp |lt(x(t n ) + Z,t n ) - ^(2;) I > . 

ze[-L,L] n ^°° 

Proof: We argue as in the proof of Proposition 13.21 If {t n } n£ N is the sequence given by 
Proposition 14. 1[ we know that (upon extracting a subsequence) u(x(t n ) + z,t n ) converges 
in Hy oc (M.) to a limit w^z) which satisfies 

Moreover, |^oo(-2)| < £ for all z > 0, ^^(z)! < B for all z < 0, and |ti>oo(0)| = £• Arguing 
as in Remark 13. 3\ one can also show that w^z) — > as 2; — ► +00. These properties 
together imply that Cqo = c* and that Woo = h e , see hypothesis H3 in the introduction. 
□ 
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Corollary 4.3 For all c > c* and all w G H^(M), one has 

e c [w] = [ e cx (~ w\xf + F(w(x))) dx > . 

Proof: Assume first that w G H^(R) for all d > 0, and that w — 1 E H\R_). If u(x,t) 
is the solution of (pQ) with initial data u(x, 0) = we know from Proposition 13.21 and 

Corollary 14.21 that the invasion point x(t) defined by ( JT7I) satisfies x(t)/t — >■ c* as t — > oo. 
Thus, for any c > c*, the quantity y c (t) = x(t) — ct converges to — oo as t — > oo, so that 
-E'c(^) > for all t > 0. In particular, £7 C (0) = £ c [itf] > 0. Letting c — > c*, we also obtain 

4>] >o. 

Assume now that c > c* and that 10 G H^(M). For any n > 1 we define 

w n (x) = w(x)x(x—n) + (1 — w (x))x(x+n+l) , x G R , 

where x G C°°(R), x(a;) = 1 for a: < and x(^) = for a; > 1. Clearly u>„(x) = 
for x G [— n, n], whereas w n (x) = for x > n + 1 and u>„(x) = 1 for x < — n — 1. Thus 
w n G i2£(R) for all c' > and w n — 1 G H\R^), so that £ c [w n ] > for all n G N by 
the preceding argument. Moreover it is straightforward to verify that £ c [u> n ] — > £ c [w] as 
n — > oo, hence £ c [iy] > 0. □ 

Equipped with these results, we are now able to prove that the solution u(x, t) con- 
verges for all times towards a travelling wave, locally in space around the invasion point. 



Proposition 4.4 For all L > we have 



L 

e c * z (u t + c*u x ) 2 (x(t) + z, t) dz ► . (47) 



Proof: We argue by contradiction and assume that there exist L > 0, 5 > 0, and a 
sequence t n — > oo such that 

J e c * z (u t + c*u x ) 2 (x(t n ) + z, t n ) dz > 5 , (48) 

for all ?i6N. Let y{t) = x(t) — c*i. If the sequence {y(t n )} n< z^ has a subsequence that is 
bounded from below, then we easily get a contradiction as in the proof of Proposition 14.11 
(case 1). So it remains to consider the case where y(t n ) — ► — oo, which requires a new 
argument. Without loss of generality, we can suppose that t n+1 > t n +T for all G N, where 
T > is as in Lemma [2T2l and that y(t n ) < —n — 1. Upon extracting a subsequence, we 
can also assume that u(x(t n ) + z,t n ) converges in i?^ c (R) towards a limit w^z). 



Given some (large) n G N, we take c n < c* such that y(t n ) = (c n — c*)t n , see Fig. 4. 
Since c n — > c* as n — > oo, we can assume that c„ > c*/2. Let y c „(£) = + (c* — c n )t, 
so that y Cn {tn) = 0. For each k — 0, 1, . . . , n we have by ([3] 



< ^(4 + T ) < e -7ib^(t fc ) + ^2 e «-(fc»(**)+Q.) 
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\C n C*)t 



y{t 3 



y(t n ) 



Fig. 4: If y(t n ) — > — oo a contradiction is obtained by considering the dissipation of the energy E Cn in a 
moving frame with speed c n < c* on the time interval [0, t n ], where y(t n ) = (c„ — c*)i n . 



hence 



(InTn p CnC ° J^- 



3=1 



(49) 



We now define fc(n) as the largest integer e N such that 

(c* - c n )tj < 1 + - , for all j = 0, 1, . . . , k . 

Since c n —> c*, it is clear that fc(n) — > oo as n — ► oo. Moreover, fc(n) < n as (c* — c n )t n 
—y{t n ) > n + 1 by assumption. For k = k(n) and j < k, we have 

y~cn(tk-j) = y(t k -j) + (c* - c n )t fc _j < —{k — j)/2 , 

hence it follows from fl49l) that 



E Cn (t n ) < E Cn (t k(n) ) < e - k ^E Cn (0) + ^^ k(n)e-^ k W-V 

c 

where 7* = min(7T , c*/4). Taking the limit n — > 00 and using the fact that £^(0) is 
uniformly bounded, we conclude that 

limsupE Cn (t n ) < . 

n^oo 

Now, since y Cn (t n ) = by our choice of c n , we have 

E Cn (t n ) = [ e c - z ( l -u 2 x + F{u)){x{t n ) + z,t n )dz , 
Jr v ^ 7 

hence taking the limit n — > 00 and using Fatou's lemma we obtain 

4>oo] = / e^dwUz) 2 + F( Woo (z))) dz < liminfE Cn (t„) < . 
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In particular G if^(R), hence it follows from Corollary 14.31 that £ Ct [woo] = 0. On the 
other hand, in view of (148p . we have 



> 



^cj^oo] = J e c * z (w'^(z) + c^w'^z) - F'(woo(^))) dz 

e c * 2 (<(z) + c^(z) - F'( Woo (z))) 2 dz 

= lim / e c * z (u t + c*u x ) 2 (x(t n ) + z,t n )dz > 5 . 

Thus, if u(x,t) is the solution of ([T|) with initial data u(x,0) = u^x), then E C:t (0) = 
£c[woo) = and E' Ct (0) = -P c Jwoo] < -5, hence £ c ,(i) < for all t > 0. This 
contradicts the conclusion of Corollary 14.31 □ 

It is now a straightforward task to conclude the proof of Proposition 11.21 Using 
Proposition 14.41 and proceeding as in Corollary 14.21 we see that u(x(t) + z,t) converges to 
Woo^z) = h £ (z) in H 2 ([— L,L}) for any L > 0. On the other hand, arguing as in (I3"9"j) . we 
find for any c G (0, c*): 

limsup / e cz (ul + u 2 )(x(t) + z,t) dz < — < oo . 

Moo Jo CK' 

This implies in particular that u(x(t) + z, t) converges to zero as z — > +oo uniformly in 
t > 0, hence u(x(t) + z,t) converges as t — > oo to h E (z) uniformly for all z G [-L, +oo). 
This proves (fl8l) . 

It remains to verify that the map t i— ► is C 1 for large i and satisfies — > c* as 
i — > oo. Using (fTSl . (120]) . and an interpolation argument, we find for any L > 0: 

sup \u x (x(t) + z,t) — h' £ (z)\ > . 



ze[-L,L] 



t—too 



As h' £ (0) < 0, this implies in particular that u x (x(t),t) is bounded away from zero for 
t sufficiently large. Since u(x(t),t) = e for t large, the Implicit Function Theorem then 
asserts that x(t) is different iable with 

x'(t) = — ^4~7~v ~\ ) f° r au sufficiently large t > . 

u x (x(t),t) 

On the other hand sup| z | <L + c#Wa.(s(t)+/Z,t)| — > as t — ► oo by ( 1201) and 

(|4"T1) . hence x'(t) — > c* as t — > oo. □ 



5 Repair behind the front 

This final section is devoted to the proof of Corollary [L3J We follow closely the arguments 
given in [131 Section 9.6], with a few simplifications. 

Let u(x,t) be a solution of with initial data uq satisfying the assumptions of 
Proposition 11.21 According to (ITS]) , we can find a time sequence t n — > oo such that 
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t n +i > tn + n + 1 for all n G N, and 

sup \u(x(t) + z,t) — ft. e (^)| < , for all t > t n . (50) 

ze[-2n,+oo) 71 + 1 

Let # : K — > [0, 1] be a smooth, nondecreasing function satisfying 8(x) = for x < and 
6{x) = 1 for x > 1. We define a map x : [0, +oo) — >• R by imposing, for all nGN: 



= - n - 0(— — -^—) , for all i G £ n+1 ] . 



It is clear that — n — 1 < < x(t) — n for all t G [t„, t n +i]- Moreover, there exists 
T > such that is differentiable for t > T, with < x'(t) < c* + 1 for all t > T. 
For later use we observe that, for any L > 0, 

sup (\u(x(t) + z,t)-l\ + \u x (x(t) + z,t)\) >0. (51) 

Indeed, since x(t) w x(t) — n for t G [i n ,t n+ i], the estimate on |w — 1| is a consequence of 
fl50|) and of the fact that ft. e (^) — > 1 as z —>■ — oo. The result for Itij.] then follows from the 
a priori bound (120]) by interpolation. 

We next consider the truncated energy function 



,2 

where F(u) = F(u) - F(l) > and 



E{t) = I <f>{x,t)(^u 2 x (x,t) + F(u(x,t))} dx 



1 if x < x(t) , 

e x(t)-x if x > x(t) . 



Since u(-,t) — 1 G i/ 1 (lR_) and t) G if 1 (IR + ), it is clear that E(t) is well-defined and 
finite for all t > 0. Moreover, i£(t) is differentiable for t > T and a direct calculation 
shows that 

£?'(*) = - f (j)(x,t)u 2 t (x,t) dx + [ <p(x,i)\x'{t)0-u 2 x + F(u)) + u x u t \ dx 
Jr Jx(t) L \2 / J 

< -1 [ (f)(x,t)u 2 t (x,t)dx + (c*+l) / e" 2 (^ + F( M ))(x(t) + z,t)d,2 . 

In view of (1201) . (I5TI) . the last integral in the right-hand side converges to zero as t — > oo. 
Since E(t) > for all t > 0, it follows that there exists a time sequence — ■> oo such that 

>(x,4K 2 (x,<Jdx > . (52) 

Now, we claim that 

sup \u(x,t' n ) — 1| ► . (53) 

xe(-oo,x(t' n )} ' n ^°° 
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Indeed, if this is not the case, there exist a positive constant e' , a subsequence {t"} n eN 
of {t^} nS N, and a sequence {x n } ne N C K such that x n < x{t'^) and |w(x n ,t") — 1| = e' 
for all )i6N. Without loss of generality, we can assume that e' > is sufficiently small 
so that the only bounded solution of the differential equation w xx — F'(w) = with 
\w(0) — 1| < ef is w = 1, see hypothesis H4 in the introduction. In view of fl5~Tj) . it is clear 
that x n — — > — oo as n — > oo. On the other hand, upon extracting a subsequence, 
we can assume that, for all L > 0, 

u(a; n + 2,0 ► ^oo(^) in H 2 ([-L,L\) , 

t^(a; n + 2,tj|) ► Woo(^) in L 2 ([-L,L]), 

n— >oo 

where G # 2 (K) H L°°(R) and ^ G L°°(M) satisfy w^z) = w'^{z) - F'( Woo {z)). 
However, it follows from (J52"j) that riioo = 0, hence Woo is a bounded solution of the 
differential equation w'^ — F'{w 00 ) = which satisfies |woo(0) — 1| = e' . This contradicts 
the assumption above on e', hence fl53|) must hold. 

Finally, if we combine (150!) and ([35)1 . we obtain 

sup + z, — /i e (2r)| > 0. 

In other words, the solution u(x, t) approaches uniformly onla translate of the travelling 
wave h £ for a sequence of times t' n — ► 00. On the other hand, the classical results of 
Sattinger [17] show that, if assumptions (j3J), (jlj) are satisfied, the travelling wave h is 
asymptotically stable with shift in the space L°°(R). In other words, Eq.(jS]) holds for any 
solution of ([T]) which is sufficiently close (uniformly on M) to a translate of h. This is the 
case for u(-,t' n ) if n is sufficiently large, hence Corollary 11.31 is proved. □ 
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